Mon. Not. R. Astron. Soc. 000, [HM] (2010) Printed 20 January 2013 (MN I*TeX style file v2.2) 



Global asteroseismic properties of solar-like oscillations 
observed by Kepler: A comparison of complementary 

- - - analysis methods 

^ ! 

(N G. A. Verner^'^*, Y. Elsworth^, W. J. Chaplin^, T. L. Campante^''^, E. Corsaro^, 

P. Gaulme^ S. Hekker^'^, D. Huber^ C. Karf, S. Mathur^, B. Mosser^o, 
5 T. Appourchaux^ J. Ballot^^'^^ T. R. Bedding^ A. Bonanno^ A-M. Broomhall^ 
R. A. GarciV^, R. Handberg^, R. New^^, D. Stello^ C. Regulol5'l^ 
I. W. Roxburgh^, D. Salabertl5'l^ T. R. White^ D. A. Caldwell^^ 
p^' J. L. Christiansen^^, M. N. Fanelli^^ 

^Astronomy Unit, Queen Mary University of London, UK 



Oh! 
I 



^School of Physics and Astronomy, University of Birmingham, UK 

^ Centro de Astroftsica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 
"^Department of Physics and Astronomy, Aarhus University, DK-8000 Aarhus C, Denmark 
^ ^ INAF Osservatorio Astrofisico di Catania, Italy 
5_j ^Institut d'Astrophysique Spatiale, UMR 8617, Universite Paris Sud, 91405 Orsay Cedex, France 

, Astronomical Institute, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands 
' ^ Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, Australia 
I I High Altitude Observatory, Boulder, CO, USA 
^°LESIA, CNRS, Observatoire de Paris, France 
, CNRS, Institut de Recherche en Astrophysique et Planetologie, 14 Avenue Edouard Belin, 31400 Toulouse, France 
^ . ''-^Universite de Toulouse, UPS-OMP, IRAP, Toulouse, France 
^ ■ Laboratoire AIM, CEA/DSM-CNRS-U, Paris Diderot, IRFU/SAp, Gif-sur-Yvette, France 

'-^Materials Engineering Research Institute, Sheffield Hallam University, UK 
l/^ , '-^ Instituto de Astroftsica de Canarias, La Laguna, Tenerife, Spain 

■ Universidad de La Laguna, La Laguna, Tenerife, Spain 
^ ■ ^"^ SETI Institute/NASA Ames Research Center, Moffett Field, CA 94035, USA 

^ Bay Area Environmental Research Inst./NASA Ames Research Center, Moffett Field, CA 94035, USA 

20 January 2013 



X 



_ ABSTRACT 

. We present the asteroseismic analysis of 1948 F-, G- and K-type main-sequence and 

subgiant stars observed by the NASA Kepler Mission. We detect and characterise 
solar-hke oscillations in 642 of these stars. This represents the largest cohort of main- 
sequence and subgiant solar-like oscillators observed to date. The photometric obser- 
vations are analysed using the methods developed by nine independent research teams. 
The results are combined to validate the determined global asteroseismic parameters 
and calculate the relative precision by which the parameters can be obtained. We 
correlate the relative number of detected solar-like oscillators with stellar parameters 
from the Kepler Input Catalog and find a deficiency for stars with effective tempera- 
tures in the range 5300 < Tcff < 5700 K and a drop-off in detected oscillations in stars 
approaching the red edge of the classical instability strip. We compare the power-law 
relationships between the frequency of peak power, t'max, the mean large frequency 
separation, Av, and the maximum mode amplitude, Amax, and show that there are 
significant method-dependent differences in the results obtained. This illustrates the 
need for multiple complementary analysis methods to be used to assess the robustness 
and reproducibility of results derived from global asteroseismic parameters. 

Key words: stars: fundamental parameters - stars: interiors - stars: oscillations - 
stars: solar-type. 
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1 INTRODUCTION 

Many stars are observed to exhibit 'solar-like' oscillations 
that are stochastically excited by near-surface convectio n 
fe.o. lBrown fc Gillilandlll994l : IChristensen-Dalsgaardll2004 ). 
These global oscillations can be observed either in spectro- 
scopic Doppler velocity measurements or by detecting small 
variations in brightness. Similar global oscillations taking 
place on the Sun have been studied in great detail and 
have provided a wealth of i nformation on the solar interior 
structure (see I Chaplin fc Basu 2008i , for a review) . Solar- 
like oscillations have been detected on a number of main- 
sequence stars using multi -site ground-based spectroscopy 
iBedding fc KieldsenI [2OO7I ) and the recent high-precision 
photometr ic observations of the French-led CO-ROT space 
telescope (jBadin et all 120071 ) . This has led to the first de- 
tailed asteroseis mic studies of solar-like osc il lations in main- 
sequence stars llAppou rchaux e t al. 200S |_ Benomar et al.l 
I2OO9I : iDeheuvels et al.l l'2010: Mat hur et ahlboiOal l. 

The NASA Kepler space telescope simultaneously mon- 
itors the brightness of over 150,000 stars with the pri- 
mary mission of detecting the characte ristic dips in bright- 
ness caused by t ransiting exoplanets (jBorucki et al.l l20ld : 
iKoch et al.ll2010l ). The high precision available in the data 
enables them to also be used for asteroseismic studies 
l|Gilliland et al.|[2010al ). The majority of stars observed by 
Kepler have photometric measurements taken at a mean 
'long' cadence of 29.4 minutes. However, a subset of up 
to 512 stars at any given time are observed on a mean 
'short' cadence of 58.85 s econdfl These high- precision, 
short-cadence observations ijGilliland et al.ir2010b[ ) are ideal 
for detecting solar-like oscillations, which have typical peri- 
ods of the order of a few minutes. 

In order to fully characterise a star using asteroseis- 
mology, it is desirable to have accurate estimates of in- 
dividual p-raode parameters - including the frequencies, 
amplitudes and lifetimes of a large number of modes for 
which the angular degree, I, and radial order, n, have 
been identified. This can be achieved using a variety 
of classical maximum-likelihoo d peak-bagging technique s 
l|Toutain fc Appourchauxl Il994l : lAppourchaux et all 119981 ) 
or si milar alternative m ethods such as a Bayesian approach 
(e.Q. lBrewer et al.ll2007l ). However, this process is often time 
intensive and only appropriate for data above a certain 
signal-to-noise level. The average asteroseismic parameters, 
indicative of the global structure, are readily obtainable us- 
ing automated analysis methods and can incorporate data 
with a lower signal-to-noise ratio for which a full peak- 
bagging analysis is not always possible. 

In this paper, we compare the global asteroseismic pa- 
rameters obtained from the analysis of 1948 F-, G- and K- 
type stars observed during the first seven months of Ke- 
pler operations and carried out using the methods of nine 
teams. This incorporates all short-cadence observations of 
stars identified as potentially showing solar-like oscillations, 
which includes a number of giants and subgiants in addition 
to stars on the main sequence. We use simulated data to 
test the precision and accuracy of each method and develop 



^ The actual cadence varies between 58.74 and 59.00 seconds once 
the observations are corrected for movement around the solar- 
system barycentre. 



a process to detect spurious outliers and ensure that data 
from different teams are self-consistent within the stated 
confidence bounds. The fraction of stars for which we can 
detect oscillations depends on the stellar parameters them- 
selves and we correlate our detections with th e parameters 
available from the Kepler Input Catalog fKIC. lBrown et ahl 
l201lh . We then derive power-law relationships between the 
verified global parameters and draw conclusions from our 
results. 



2 DATA 

We have used the short-cadence data from the first seven 
months of observations made by Kepler. These data consist 
of 2229 light curves obtained from stars of Kepler-haxiA mag- 
nitude Kp = 6.9 to 12.1, each of approximately one-month 
duration and obtained with a mean cadence of 58.85 seconds. 
These stars were predetermined to be potential solar-like os- 
cillators based on their characteristics in the KIC. Some of 
the stars were observed during multiple one-month runs, 
resulting in 1948 unique stars in the complete set. This rep- 
resents 74 % of the stars observed at short cadence during 
the initial survey phase of the Kepler Mission. 

The li ght curves we r e first prepared in the manner de- 
scribed in iGarcia et all (|201ll ) starting with the raw data 
provi ded from the Kepler Mission pipeline (jjenkins et al.l 
l2010h . By using a common pre-processing method to remove 
spurious points and detrend each light curve, we ensured 
that all of the analysis teams started with the same input 
data for their methods. We had found that variations in 
the treatment of outlying data points and different meth- 
ods used to interpolate missing points or remove long-term 
trends introduced a significant bias into the asteroseismic 
parameters obtained for each star. 



2.1 Global parameters 

The prepared time series and corresponding power spectra 
were analysed to give estimates of the asteroseismic param- 
eters sensitive to the global structure of the observed star. 
The most prominent feature of an acoustic power spectrum 
is the near-regular spacing between modes of the same an- 
gular degree and consecutive orders. The average of this so- 
called large frequency separation {/S.v) gives an estimate of 
the total acoustic radius of the star, i.e. the sound travel 
time from the centre to the surface: 
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The average large frequency separation is also known to scale 
with the square root of the mean density of the star (e.g. 
IChristensen-Dalsgaard|[l993l ) . giving the scaling relation 
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The heights of the p-mode peaks in a power spectrum 
are modulated by an envelope that attains its maximum 
value at a frequency we denote as Vma.x. It is usual to 
make the assu mption that z^max scales with the acoustic cut- 
off frequency l|Brown et al.lll99ll : lKieldsen fc Beddindll995l : 
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1 Mathur et al. (2010b') : ^ Mosser fc Appourchaux (2009) : ^ Hekker et al. (2010); Bonanno et al. (2008) : 
5 Verner fc Roxburgh (2010) : ^ Campante e t al. (2010); Karoflt et al. (2010) : ^ Huber et al. (2009) : 
t only used for the analysis of simulated data. 



Table 1. Brief description of the methods used by each analysis team. Abbreviations used are: SPS - smoothed power spectrum, PSPS 
— power spectrum of the power spectrum, TSACF - ti me series autocorrelation function, PSACF - power spectrum autocorrelation 
function, K08 - method based on lKieldsen et al.l l|2008lJ) . 



iMosser et aLll2010[ ). which when determined for an isother- 
mal stellar atmosphere gives a scaling relation in terms of 
mass, radius and effective temperature of 
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The height attained by a single p-mode profile is related 
to the intrinsic amplitude of the underlying stellar oscilla- 
tion but is also affected by the lifetime of the mode. Modes 
with shorter lifetimes have profiles with greater linewidths, 
spreading power over many frequency bins. If the power 
spectrum is converted to power spectral density by dividing 
by the frequency resolution, the intrinsic RMS amplitude of 
a radial mode can be obtained as 
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where Hg is the height of the underlying limit-spectrum pro- 
file, usually modelled by a Lorentzian function, and F is the 
corresponding linewidth. To obtain the intrinsic amplitude 
of non-radial modes, an extra factor appears in Eq.Uto ac- 
count for the averaging of the oscillation eigenfunction over 
the visible stellar disk and the effects of limb darkening. The 
amplitude that would be attained by a radial mode with a 
frequency of Vmax we refer to as Amax and reflects the ex- 
citation mechanism present in the star. If we assume an 
adiabatic stellar atmosphere, the value of j4max, once cor- 
rected to a bolometric amplitude, is expected to scale with 
the RMS velocity perturbation (uosc) a nd inversely with the 
square root of effective temperature (|Kieldsen fc Beddind 
1 19951 ): 
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The surface layers of a star introduce a significant 
frequency-dependent per turbation into the m easured fre- 
quencies of p modes (e.g. iKieldsen et al.ll2008al ). The modes 
of degree I = and I = 2 are closely spaced in frequency 
and are therefore affected similarly by the surface layers. 



The small frequency separation between these modes ((5 1^02) 
is sensitive to the physical conditions in the deep interior 
and is therefore an important indicator of stellar structure, 
age and evolutionary state. Additionally, the ratio of large to 
small frequency separations has been shown to discriminate 
well between stellar evolutionary models and is insensitive 
to variations in the surface layers l|Roxburgh fc Vorontsovl 
I2OO3I : iRoxburghllioOsI , I2OI0I ) . 



3 VALIDATION AND VERIFICATION OF THE 
DIFFERENT METHODS 

The prepared time series were analysed by the automated 
methods of nine teams. Each method returned some or all of 
the globally averaged seismic parameters outlined above for 
each star. Although the exact methods used by each team 
varied, many of the methods used s imilar techniques based 
on autocorrelation a nalyses (e.g. [Roxburgh fc Vorontsovl 
I2OO6I : lRoxburghll2009l ) to obtain Av and i^max o r the po wer 
spectrum smoothing method of IKieldsen et ahl (|2008lJ ) to 
obtain Amax and Vmmx- The methods used are briefly out- 
Uned in Table [T] 



3.1 Simulations 

In order to test the performance of a representative sample 
of the analysis methods, a set of 309 simulated Kepler-like 
time series each of one-month duration were prepared and 
analysed by six teams (A2Z, COR, KAB, OCT, QML, SYD) 
in a blind hare-and-hounds exercise. The stellar parameters 
were chosen to cover a range in values of Ai/, fmax, ^max 
and Svq2 in accordance with the expected parameters of the 
stars showing solar-like oscillations in the Kepler field. The 
values returned by the teams were compared with the in- 
put values of asteroseismic parameters for each simulated 
star to determine the bias, precision, number of outliers and 
accuracy of formal uncertainties for each parameter. Due 
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Figure 1. Histograms showing the relative deviation of the fitted 
values of Ai/, z^max, Amax and from their true values as 

determined from the simulated data results across all methods. 

to the large variation of the parameters over the stars in 
the cohort, and the assumption that the inverse granula- 
tion timescale and wi dth of the p-mode power excess both 
scale with frequency (|Huber et al.ll200 9|), we have assumed 
that the formal uncertainties scale with frequency also and 
therefore used relative uncertainties in our analysis. The his- 
tograms showing the relative deviations of each parameter 
from the corresponding input values across all methods are 
shown in Fig. m 

The simulated data were generated using realistic 
stellar-model frequencies as input parameters, thus the val- 
ues of A:/ and 5vo2 for each simulated star vary with fre- 
quency. The comparison reference values of Aj/ were deter- 
mined using a cubic spline interpolation to the spacing of 
the seven radial modes in the vicinity of t'max as a function 
of frequency, evaluated at i^max- The reference values for 5vo2 
were obtained in a similar way. However, the (5i'o2 spacing, 
particularly for evolved solar-like oscillators, may not vary 
smoothly with frequency due to mode 'bumping' caused by 
coupling with g modes confined within the stellar interiors 
(White et al., in preparation). Stars in the simulated set 
which did not have a smoothly varying 5i^o2 spacing in the 
region around Vmax were not included in the 5uo2 comparison 
analysis. 

3.1.1 I^v 

The analysis of simulated data showed that Ai/ was the 
most precisely determined parameter. For the results of each 
method, we identified outliers as those points lying outside 
of a relative range of 5 a from the true value. We used the 
median absolute deviation, scaled by 1.4826 assuming a nor- 
mal distribution, for our estimate of ct to reduce the influ- 
ence of outliers on our measure of the spread. The average 
relative number of outliers we found across all methods was 
6%, although this varied from 1 % to 11%. With the out- 
liers removed from the results of each method, the average 
standard deviation of the remaining distributions was 1.2 % 



and varied between 0.7% and 2.0%. We applied the t-test 
to the distributions and found that four (out of twelve) had 
centroids sufficiently far from zero that the null hypothesis 
(mean is zero) could be rejected at the 95 % level, hence in- 
dicating a small bias. The centroids of these distributions 
were all between 0.3% and 0.4%, indicating that any over- 
estimation bias is at a level that is well within the expected 
1 (7 error distribution. These biases were not apparent in the 
later analysis of Kepler data. 



O. 1 . i^raax 

The simulations showed that the spread in values of t'max 
is wider than that for Ai/. The values for comparison were 
determined from the frequency at which the amplitude en- 
velope of radial modes attained its maximum value. In the 
simulations, as in real stars, the non-radial / — 1 modes have 
a slightly higher visibility than radial modes due to their in- 
creased spatial contribution when averaged over the visible 
stellar disk and the effects of limb darkening. The ampli- 
tude envelope for non-radial modes is assumed to peak at 
the same frequency as that for radial modes, however the 
determined value of fmax may be sensitive to which modes 
are present in the region around (--max. Many of the meth- 
ods use a wide smoothing filter in the power spectrum to 
minimise this effect. 

We identified outliers in the same way as for Af and 
found that the average proportion of outliers was 5 %, vary- 
ing from 1 % to 9 % over the methods. Once these outliers 
were removed, the average standard deviation of the re- 
maining distributions was 3.0% and varied between 2.1% 
and 4.3 %. Applying the t-test to the distributions indicated 
that three (out of eight) had centroids that rejected the zero- 
mean hypothesis at the 95 % level. The mean value of these 
distributions were —0.9%, 1.1% and 1.2%, the largest of 
which corresponded to 0.4 a. 

The presence of a bias in the estimation of fmax may 
be due to the model used by each method to fit and remove 
the granulation background or the function used to fit the 
p-mode power excess. The methods that were identified as 
having a small bias in the analysis of simulated data were 
subsequently improved upon before the analysis of real data. 

3.1.3 ^max 

The results obtained for Amax showed a much wider spread 
than for Av and i^max. The relative number of outliers was 
slightly lower than for the other parameters due to the in- 
crease in the measured 1 a spread of values, with an aver- 
age of 3 % and varying between 1 % and 6 %. The outlier- 
removed mean standard deviation over all methods was 
14.4%, and varied between 9.7% and 15.6%. 

By applying the t-test to each distribution, the cen- 
troids of only one distribution (out of eight) was sufficiently 
close to zero to be consistent with the zero-mean hypothesis 
at the 95 % level. The other distributions had mean values 
between —15 % and 2 %, which indicated an underestimation 
bias in all cases except one. Despite the presence of a bias, 
all of the methods showed a strong correlation between the 
values of Amax determined and the simulation input values, 
with correlation coefficients between 0.92 and 0.96. Fitting a 
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Figure 2. Histograms siiowing the distributions of tlie relative deviation of Amax from the correct value for each method determined from 
the analysis of simulated asteroseismic data. The method names and number of artificial stars included are given in the titles. Dotted 
lines show the distributions based on the Amax values returned by each method, solid lines show the distributions after a calibration 
scaling to match the true values. 



linear relation to the variation of the determined Amax as a 
function of the input values suggested that a simple scaling 
is sufficient to calibrate the amplitudes determined by each 
method. 

Determining the maximum mode amplitude in an auto- 
mated way requires assumptions to be made about the rela- 
tive contribution of each mode deg ree to the amplitude p er 
radial order {e.g. the c parameter in lKieldsen et al.l(|2008bl )l. 
the amount of smoothing applied to the power spectrum 
(a wider smoothing window can reduce the amplitude of 
the p-mode power excess) and also an appr opriate contri- 
bution to model t he stellar background (e.g. lHarvevlll985l : 
iKaroff et al.ll201ol ). As each method addressed these prob- 
lems in a slightly different way, the amplitudes they obtained 
were scaled differently from each other. To calculate a cali- 
bration factor for each method, a linear regression was per- 
formed on the Amax values from each team against the true 
values for the simulated data. The Amax values were then 
scaled by this calibration factor before any further analysis. 
The relative deviations of the Amax results from each method 
before and after the calibration was applied are shown in 
Fig.H 



3.1.i 5vo2 

The values of the Svm small separation were only determined 
by two of the methods (COR, QML). Repeating the same 
outlier detection procedure applied to the other parameters 
indicated that the relative number of outliers was 17 % and 
8 % for these two methods. The average 1 a spread in values 
of 5vo2 was the largest of all the parameters at 26 %, reflect- 
ing the difficulty with which one can precisely determine this 
small spacing in an automated way from data of just one- 
month duration. Despite the wide distributions, applying a 
t-test found that both methods had distribution centroids 



that were consistent with the zero-mean hypothesis at the 
95 % level, indicating that no bias was present. 

3.1.5 Formal uncertainties 

The formal uncertainties of each parameter from each 
method were tested for consistency against the known devi- 
ations from true values using the test. In many cases the 
formal uncertainties were found to be inconsistent with the 
accuracy of the results. To obtain a per degree of freedom 
of unity for the parameters from each method would require 
scaling the formal uncertainties by a significant amount. 
Despite this inconsistency, the relative formal uncertainties 
were correlated with the relative absolute deviations of the 
results from their true values, particularly in the case of /S.u 
(with a correlation coefficient of 0.30) and t'max (0.29). The 
formal uncertainties and relative deviations of Amax were 
less strongly correlated (0.20), however all of these correla- 
tions are sufficient to reject the null hypothesis (i.e. random 
points with zero correlation) at the 0.1 % level. In the case 
of Svm , the stated uncertainties and deviations were not sig- 
nificantly correlated, however the sample size was relatively 
small because only two of the methods returned results on 
5uo2 for the simulated stars. 

3.2 Application to Kepler data 

The prepared Kepler Hght curves were analysed using the 
methods of all nine teams (see Table [l}. The methods did 
not all detect oscillations on the same targets so there was 
some non-overlap in the sets of stars for which oscillations 
were detected by each method. In order to maximise the 
number of stars for which we obtained verified global aster- 
oseismic parameters, it was desirable to combine the results 
from all methods rather than ranking them in some way 
and selecting only the best. The methods used previously 
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Figure 3. Histograms showing the distributions of the relative deviation of Amax from that given by the SYD-1 method for Kepler 
data. The method names and number of stars included are given in the titles. Dotted lines show the distributions before calibration by 
peak-bagging results, solid lines show the distributions using the corrected amplitudes (including those by SYD-1). 



to analyse the results from simulated data required knowing 
the correct parameters a priori. In the case of real data we 
constructed alternative methods to detect outliers, calibrate 
the maximum mode amplitudes and test the consistency of 
formal errors. 

3.2.1 Outlier rejection 

The results from the simulated data showed that for all of 
the parameters, there were significantly more outliers at 5 a 
than might be expected assuming a normal error distribu- 
tion. This is an indication that the error distribution is not 
solely random and has a dependence on the varying signal- 
to-noise ratio of the stars present. It also highlights the lim- 
itations of the automated analysis methods in cases where 
the signal-to-noise ratio of the data is low. 

As a verification procedure, we required that for each 
parameter determined for each star the results from at least 
two independent teams were contained within a range of 
fixed relative size centred on the median value. Results out- 
side of this range were iteratively removed until either all re- 
sults were in agreement or fewer than two results remained. 
The range needed to be small enough to reject outliers but 
not so small that results were removed due to the inherent 
scatter in the data. 

The range size for each parameter was defined using 
the results from the simulated data. We used the smallest 
range such that the outliers identified at the 5cr levels de- 
scribed previously would also be detected using the fixed- 
range method. The relative ranges that we calculated were 
±3.5 % (Ai/), ±10.5 % (z^max) and ±21.5 % (yl^ax and 5vo2). 

3.2.2 Amax calibration 

To identify any significant biases in the asteroseismic param- 
eters derived from Kepler data, we compared the results of 
each method to those of a reference method (SYD-1). As was 



the case for the simulations, the parameter with problematic 
biases was Amax. In order to calibrate the amplitudes of each 
method, the maximum mode amplitudes of 22 high signal- 
to-noise stars within the cohort were determined using an 
established peak-bagging technique and used as a reference 
set to obtain Amax correction factors for each method. The 
peak-bagging technique was tested using a subset of 48 sim- 
ulated stars with high signal-to-noise ratios and found to 
give an unbiased estimate of Amax to a relative precision 
of 6.2%. The relative differences between the ylmax values 
for Kepler data determined by each method and that given 
by the reference method before and after calibration by the 
peak-bagging results are shown in Fig. [S] This illustrates the 
improved agreement in the values determined for ylmax by 
each method after calibration. 



3.2.3 Formal uncertainty correction 

In the case of real data, it was not possible to calculate 
uncertainty-scaling factors using the test without know- 
ing the correct values a priori and so a method was devised 
to estimate these using the scatter in the results obtained for 
the same stars from multiple teams. For those stars where at 
least five teams gave results that were in agreement within 
the thresholds defined above, the standard deviation of these 
results was used as an estimate of the expected uncertainty. 
A scaling factor was then determined from the ratio of this 
standard deviation to the formal uncertainties stated by 
each team. Single average scaling factors for each method 
were then calculated as the median of the individual scaling 
factors. The distributions of the uncertainty-scaling factors 
are shown in Fig. |4] and the median values for each method 
are given in Table [2] The distributions in Fig. |4] are encour- 
aging and show that once scaled by a single scaling factor, 
the uncertainties of each method agree well with the spread 
observed between results obtained from multiple methods. 
The error correction procedure ensured that the data 
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Figure 4. Distributions of uncertainty-correction factors {i.e. standard deviation divided by formal uncertainty) for Au, Vmax and 
j4-max- The correction factors for each method have been scaled by their median values (Tabled and combined in the histograms for 
each parameter. The error bars show the standard deviations in each bin when determined using individual methods. 



were self-consistent within stated uncertainties, however 
there is an additional internal uncertainty which must be in- 
cluded. For a given star we expect the scatter in the results 
obtained from multiple teams to be less than a true measure 
of the uncertainty because all of the methods were analysing 
the same data. If we use the scatter by itself to calculate 
the correction factors then they will be underestimated. By 
estimating the uncertainty-scaling factors using the scatter 
method on the simulated data and then applying the test 
to check the consistency of the scaled uncertainties, we de- 
termined the additional internal uncertainty-scaling factors 
required for Av (1.20), i^max (1-19) and Anax (l-H)- 

For certain sets of results, the scaling that must be ap- 
plied to the formal errors to ensure they are consistent with 
multiple alternative methods is considerable (see Table [2]). 
For cases where this scaling is large {i.e. greater than four), 
the reliability of the scaled uncertainties is questionable. In 
addition to this, some of the methods did not calculate any 
formal uncertainties on their results. For these data, a con- 
stant relative uncertainty was assumed to be more appro- 
priate than using the rescaled values or a constant absolute 
uncertainty. The size of these indicative uncertainties was 
determined from the median relative standard deviation of 
the results of each parameter where at least five teams gave 
consistent results and the additional internal uncertainty 
contribution described above. This gave an indicative rel- 
ative uncertainty of 1.8% for Av, 3.8% for t'max and 9.8% 

for j4max. 

Where these results were used in subsequent analysis to 
derive model coefficients, we repeated the calculations with 
both scaled and indicative errors and found that the coef- 
ficients derived were in agreement at the level of precision 
available in the data. 



4 DISTRIBUTION OF PARAMETERS AND 
DETECTABILITY 

The ylmax correction, verification procedure and formal un- 
certainty scaling was applied to the results of the Kepler data. 
obtained from nine independent analysis methods. The re- 
sults from stars that were observed in multiple one-month 
runs were averaged to give a single set of results for each 
star. Of the 1948 unique stars in the cohort, 642 (33 %) were 



Method 


(TAi^ factor 


"'I'ma.x factor 


""Amax fS'CtOr 


A2Z 


0.32 


28.5 


1.12 


AAU 


0.94 


0.47 


0.75 


COR 


1.14 


3.18 




IAS 


7.55 


305.0 


* 


KAB-1 


0.30 


6.92 


8.85 


KAB-2 




2.48 




OCT-1 


2.03 


1.09 


0.74 


OCT-2 






1.07 


OCT-3 






0.89 


ORK 




* 




QML-1 


0.64 


0.89 


0.22 


QML-2 


0.45 






SYD-1 


0.93 


0.86 


1.18 


SYD-2 




0.77 


0.88 



Table 2. Median uncertainty-correction factors for each method. 
Values greater than unity indicate where formal uncertainties 
have been underestimated. Asterisks denote cases where results 
were provided without formal uncertainties. For methods where 
a scaling factor greater than four is required, we have assumed 
these uncertainties to be unreliable. 



confirmed as solar-like oscillators with their determined as- 
teroseismic parameters consistent with the results of at least 
two independent teams. The distributions of the average as- 
teroseismic parameters for each star are shown in Fig. [5] 

The stellar parameters for the stars observed by Kepler 
are provided in the KIC. This includes the apparent mag- 
nitude in the Kepler band, effective temperature, metallic- 
ity and surface gravity derived from Sloan photometry of 
pre-flight observations. The KIC effective temperatures are 
found to differ from those determined from spectroscopic 
studies which are typically 200-300K higher than those in 
the KIC. A recalibration of the KIC Sloan photometry was 
performed by Pinsonneault et al. (in preparation) who find 
effective temperatures that are ~280 K higher than the KIC 
temperatures for solar-like stars. We used these recalibrated 
temperatures with the apparent magnitude, metallicity and 
surface gravity from the KIC to compare the distribution of 
stellar parameters in the confirmed solar-like oscillators to 
those in the complete set (Fig. |6} . 

We find an expected correlation of the relative number 
of confirmed solar-like oscillators with the apparent magni- 
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4500 5000 5500 ^000 6500 7000 7500 4500 5000 5500 6000 6500 7000 7500 3 35 4 45 5 0> 

Figure 6. Left panels show the distributions of stellar parameters in the full set (open bars) and for confirmed solar-like oscillators (filled 
bars). Right panels show the ratio of these distributions with error bars calculated from Poisson statistics. 




Figure 5. Histograms of the average asteroscismic parameters 
confirmed by at least two methods. The number of parameters 
in each range is scaled by the total number to give the relative 
number in each bin. The bin widths are 5 /.tHz (A;^), 100 /iHz 
(I'max), Ippm (A 

max 

) and 0.01 (5vo2/^u). 



starting at ~6700 K with fewer than 10 % of the stars with 
Tcft > 7000 K having confirmed solar-like oscillations. This 
high-temperature drop is expected and is in agreement with 
the location of the red edge of the classical instability strip. 
We also observe a slight decreasing trend in the relative 
number of oscillators with increasing metallici ty and surface 
gravity, which is contrary to what is expected JSarnadi et al.l 
I2OO7I : iFrandsen etallbOOTl : IStello fc GiUilandll2009l ). We find 
no significant correlation between metallicity and apparent 
magnitude that would explain this trend for the stars with 
detected oscillations. However, it must be stressed that the 
KIC values for metallicity and surface gravity are assumed 
to be accurate only to ~0.5 dex. 

The average values of i^max and Ai^ are plotted against 
effective temperature in the asteroseismic H-R style dia- 
grams in Fig. [T] It is clear from these diagrams that the 
cohort includes a population of subgiants as well as main 
sequence stars, with the majority of main-sequence stars be- 
ing hotter than the Sun and with lower values of i^max and 
Af. The stars with asteroseismic parameters closer to those 
of the Sun also tend to be the brighter stars in the set. This 
is illustrated in Fig. O which shows the median values of 
Ai/, fmax and j4max as a function of apparent magnitude. 



tude of the star. At Kp ~ 9.5, we find approximately 50 % of 
the stars to be confirmed solar-like oscillators compared with 
only 20 % at Kp ~ 11.5. The effective temperature distribu- 
tions of the complete cohort and that of the confirmed solar- 
like oscillators are bimodal, corresponding to a population of 
subgiants between ~5000 and ~5500 K and a hotter group 
of main-sequence stars. There is a notable reduction in the 
proportion of confirmed solar-like oscillators in the region 
between these distributions, for 5300 < T^a < 5700 K. This 



has also been noted by I Chaplin et al. (|20l3) who suggest 
that this absence may be due to evolutionary effects that 
cause an increase in surface magnetic activity and there- 
fore reduces the detectability of oscillations. The effective 
temperature dependence also indicates a high-temperature 
fall-off in the proportion of confirmed solar-like oscillators 



5 RELATIONSHIPS BETWEEN GLOBAL 
PARAMETERS 

5.1 i-'max and At* 

It has been observed that Av and i^max follow a close 
power-law relations hip both for main sequence stars (e.g . 



Stello et al. 2009a ) and red giants (iHekker et al 



Bedding et al."2010|; iHuber et all l2010l : iMosser et al.l 



2oog 



201c 



Hekker et al. 201L). By combini ng the scaling relations Eq.[3] 
and Eq. [2l IStello et al.l l|2009ah showed that the expected 
relationship between Av and i/max is almost independent of 
luminosity and only weakly dependent on mass and effective 
temperature. To a good approximation, we can assume 
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Figure 7. Asteroseismic population diagrams produced using the average values of Au {left) and i^max {right). The subgiants are grouped 
in the upper right of the diagrams, with the main sequence stars grouped in the lower left. The diamonds indicate the position of the 
Sun. The points are coloured according to the corresponding apparent magnitude taken from the KIC. 
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Av OC (//max)" • (6) 

Using a set of 55 main sequence and giant stars, IStello et al.l 
l|2009al ) obtained a power law exponent of a = 0.772 ± 
0.0 05. This can be co mpared with the value obtained 
by iMosser et al] (|2010l '). using data from 930 red giants 
observed by the C O ROT satellite, of a = 0.75 ± 0.01. 
iMosser et al] (|2010l ) point out that the apparent discrep- 
ancy between these values relates to the different temper- 
atures and evolutionary states of the stars used in each 
study. iHekker et all (|2009l ) found a slightly higher value for 
COROT red-giant stars of a = 0.784 ± 0.003, however this 
was obtained after scaling by solar values and fitting a power 
law of the form 



which fixes the proportionality constant in Eq. E]and forces 



t he fit through the solar values of Az/ and i^mav. lStello et all 
( i2009al ) showed that an extra scaling factor appears in Eq.[7] 
that is close to unity but has a weak dependence on mass 
and TcH, therefore we expect a spread about this power law, 
as discussed in H uber et al. 

I dioio). As the set of stars used 
in this analysis cover a wide range in Tofi, and each method 
returns results on different set of stars (i.e. covering different 
ranges in Teff), we have scaled fmax and Ai^ by their solar 
values and fitted power laws of the form given in Eq. [T] As 
this restricts the power law to one degree of freedom, the ex- 
po nents obtained are n ot dir ectly comparable to t hose found 
by IStello et all (|2009al ) and IMosser et all (l2010l'). but may 
be com pared with that found for red giants bv lHekker et al.l 
(|2009t ). 

We have determined the power law exponent a from 
the results of each of the methods using a weighted least- 
squares approach. This was done using the scaled formal un- 
certainties returned by the methods and also using constant 
relative indicative uncertainties estimated from the mean 
scatter in the results and the internal uncertainties deter- 
mined from simulated data. The power law fits were per- 
formed over the full range in fmax and also individually in 
three ranges: i^max < 750 ^Hz, 750 < i^max < 1500 /xHz and 
J^max > 1500 ^Hz. The power law fits are shown in Fig. [9] 
and the derived exponents plotted in Fig. [10] and given in 
Table IXn 

Taking the weighted average of all methods when we 
used the corrected uncertainties gave a value of a = 0.795, 
with a mean error of 0.005 and a standard deviation of 0.007 
across the different methods. When we repeated this using 
fixed relative uncertainties for all of the methods we ob- 
tained a value of a = 0.794, with mean error 0.005 and a 
standard deviation of 0.008. It is reassuring that the average 
value of a is not significantly affected by the way in which 
we treat the formal uncertainties. However, in both cases 
the spread in values obtained for a by the different methods 
was greater than the average error derived for the parame- 
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Figure 10. ZXzv* vs. t^max power law exponents (Eq . Vf\ determined 
over the Vmax ranges given in eacli plot. Points with crosses are 
determined using the scaled formal uncertainties submitted by 
each group, circles are determined using constant relative uncer- 
tainties. Dotted lines shows the mean values calculated from all 
methods in each range. Methods used are given in the format 



ter. This apparent disagreement is due to the fact that each 
method gives results on a different sample of stars, sampling 
a different range of metallicities, effective temperatures and 
evolutionary states. As the degeneracy of Eq. [7]is known to 
be affected by stellar composition, it is not surprising that 
each method gave slightly different values of a. As a refer- 
ence value for the uncertainty on a we suggest that it is more 
appropriate to take this scatter into account rather than to 
use the error obtained from a single method. The i^max-Az^ 
relation is discussed further by Huber et al. (in preparation). 



5.2 



and An 



The calibrated maximum mode amplitudes from each 
method (Amax.obs) were scaled to their bolometric equivalent 
(^max = Cif_boi^max,obs) usiug the bolometric correction 
derived from the spectral response of the Kepler passband 
(|Ballotll201lh . 



CA'-boi = 1 + ai(roff — To) -f a2{Tc« - To) , (8) 

with To = 5934 K, ai = 1.349 x 10"* K"^ and a2 = -3.12 x 
10"^ K-^ The average bolometric maximum mode ampli- 
tude as a function of i/max is shown in Fig.|9l The relationship 
between Amax and i/max is expected to have a colour depen- 
dence (Eq. O , therefore power law fits were performed over 
three ranges in T^s- T^s < 5560 K, 5560 < T^s < 6250 K 
and TcS > 6250 K. The fits were performed according to 



(9) 



^max,0 



The power laws were determined using the corrected formal 
uncertainties and also the constant relative uncertainties in- 
dicative of the spread in the results from each method. The 
power law indices determined for each method are given in 
Table Ell 

Following the sca ling relations giv en in 

iKieldsen fc Beddind (| 19951 ) and ISamadi et all (|200'iD . 
for an adiabatic stellar atmosphere the maximum mode 
amplitude is expected to scale as 



An 



A 



max,0 



M/Me 



T, 



cff,0 



(10) 



where bolometric intensity variations are observed. For ob- 
servations obtained in a narrow wavelength band, the wave- 
length dependence of the measured mode amplitude results 
in a modification of the effective temperature exponent in 
Eq. [To] from —0.5 to —1.5. The assumption of an adia- 
batic stellar atmosphere is a si mple approximation and, us - 
ing narrow-band observations, iKield sen fc Bedding! (|l995l ) 
showed that an effective temperature exponent of —2.0 gave 
a better fit to t he data in the i r sarn ple. For main-sequence 
solar-like stars, ISamadi et al.l l|2007l ) calculated a value of 
s — 0.7 based on numerical models. However, using the bolo- 
metric amplitudes of red giants observe d by CO ROT and an 
effective temperature exponent of — 0.5. lM'osser et all (|201G| ) 
obtained a value of s = 0.89 ±0.02. It is apparent from these 
results that there is an additional temperature dependence 
in Eq. 1101 Here we incorporate any additional temperature 
dependence into the exponent s and fix the effective tem- 
perature exponent to —0.5. 

Rewriting luminosity in terms of radius and T^s and 
using Eq. [3] leads to: 



M/Mc, 



Tcff 



t^max, 



(11) 



The value of s can then be found by combining and rearrang- 
ing Eas. flOlandfTTIor by fitt i ng a power law to A max (Teff)"'^ 
vs. L/M (jHuber et al.lf201ol : iMosser et al.ll2010[ l. The power 
law fits and the variation in s with T^g are shown in Fig. 

We find a clear dependence of s on Tcff for all of the 
met hods. The value of s for cooler stars is in agreement 
with lMosser et all l|2010l ) who find a value of s = 0.89 ±0.02 
in a set of red giants with an average Toff of 4500 K and 
explains the different values of s obtained for main-sequence 
and giant stars. The value of s for Tcft,Q = 5777 K and the 
gradients determined from Fig. [5] are given in Table IA2I 
Taking the weighted mean values across all methods gives 



s = 0.64 - 1.18 



Tcff — Tcff .( 



T 



cfr,0 



(12) 
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Figure 9. Power law fits for stated parameters, solid lines show the fits derived from different methods. For the upper right and lower 
left panels, the fits shown are those performed over the range T^g > 6250 K. Data points shown are the mean values determined using 
all methods. 



It is not surprising that there is a remaining temperature 
dependence in Eq.[Tn]as the scahng relation rehes on a very 
simple adiabatic descr iption of the stellar a tmosp here. For 
further discussion, see iKield sen fc Beddind (|201ll ). 



6 SUMMARY AND CONCLUSIONS 

We have detected solar-like oscillations in 642 F-, G- and K- 
type main-sequence and subgiant stars in the Kepler field- 
of-view and characterised them by their large frequency sep- 
aration, frequency of peak power and maximum mode am- 
plitude. By combining results from the analysis methods of 
nine independent teams, we have verified the detections, re- 
jected outliers and devised a method to ensure the results are 
consistent within an accurate uncertainty range. It is appar- 
ent that the formal uncertainties returned from automated 
analysis methods are often inconsistent with the actual pre- 
cision of the results. Obtaining an accurate uncertainty on 
the global asteroseismic parameters is essential when using 
such results to model stellar structure. 

Using the standard deviation of the results obtained 
by multiple methods and the intrinsic uncertainty esti- 



mated from the analysis of simulated data, we have deter- 
mined the average expected relative precision of Ai/ (1.8 %), 
t^max (3.8 %) and ^m ax (9 .8 %). In the case o f red giants, 
iBedding et al] (|2010l ) and iHuber et al] (|2010l ) have shown 
that, when radial modes can be identified, the 5vq2 separa- 
tion and the e factor o f the Tassoul a symp totic law scale 
with Az^. Furthermore, iMosser et al.l (|201ll ) have demon- 
strated that the red-giant oscillation frequencies obey a uni- 
versal pattern that helps to measure A;/ very precisely and 
to identify unambiguously the mode angular degree and ra- 
dial order. For solar-like stars, the amplitudes of oscillations 
are much smaller than for red giants and the linewidths 
of the modes are often broader, reflecting th e contribution 
of increased temperature on mode damping (|Chaplin et al.l 
120091 ). We have also observed mixed modes and avoided 
crossings in evolved solar-type stars, which introduce discon- 
tinuities into the small frequency separation. These effects 
combine to make the process of accurately estimating the 
average small frequency separation a more difficult task for 
solar-like stars than for red giants. Using the analysis of sim- 
ulated data, we have shown that the average expected rela- 
tive precision oi5vo2, as obtained using automated methods, 
is only ~26%. However, for stars with a sufficient signal-to- 
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noise ratio, this can be much improved upon by using a full 
peak-bagging approach and obtaining individual p-mode fre- 
quencies. 



By characterising the stars by A.v or 



we have 



clearly identified the separate main-sequence and subgiant 
populations (Fig.[7|. The distributions of these asteroseismic 
parameters in the observed main-sequence stars show that 
we typically find smaller values of /S.u and fmax than we find 
for the Sun (i.e. Ai/q ~ 135 /xHz, t'max,© ~ 3100 ^Hz). This 
reflects the increased amplitude of oscillations in stars with 
lower i/max and the higher intrinsic luminosity of such stars 
(Fig.®. 

We have correlated the fraction of stars for which we 
detected oscillations with the stellar parameters from the 
KIC and have found a significant reduction in the propor- 
tion of solar-like oscillators with effective temperatures close 
to Teff ~ 5500 K, the temperature range that separates the 
distributions of main-sequence st ars from subgiants. T his is 
in agreement with the findings of lChaplin et al.l |20l3) who 
suggest that stars in this region may show interesting evo- 
lutionary effect s in their stellar dynamos (iGillilandl llQisI : 
iDall et al]|2010l ). which would manifest in the surface mag- 
netic activity and, therefore, decrease the detectability of 
oscillations. 

We have determined the exponents of power-law rela- 
tionships between verified A;/ and fmax values and have 
found that the scatter in the results is at a level greater than 
the formal uncertainties obtained for the individual methods 
(Fig. [To] and Table IA1[| . This illustrates the importance of 
using results from multiple methods when testing such re- 
lationships between parameters. Taking a weighted average 
of the exponents, we have determined a power law exponent 
of a = 0.795 ± 0.007, which is significantly higher than that 
found using data from a similar number of red giants using 
the same method (Hckkcr ct al. 200^), a = 0.784 ± 0.003. 
This discrepancy can be explained by the different tempera- 
tures and evolutionary states of the stars in each study and is 
also s een in the different results obtained for main-sequence 
stars fetello et al. l|2009ah and red giants l|Huber et al.ll2010l : 
[Mosser et al. 201oiv 

The maximum mode amplitudes obtained by each 
method, once calibrated by those found for stars with a high 
signal-to-noise ratio, have been used to determine the expo- 
nent s relating the dependence of mode amplitude on L/M 
(Eq.lll|l. Various authors have stated values of s between 0.7 
and 1.0 using the amplitudes calculated from main-sequence 
and giant stars. We have used bolometric-corrected ampli- 
tudes with a temperature dependence exponent of —0.5 and 
found that a strong temperature dependence remains in s. 
We have constrained this dependence by a linear relationship 
(Eq. I12|l which gives a value of s for stars in the red-giant 
regim e which agrees with that determined in iMosser et al.l 
l|2010t ) who used the same expected temperature dependence 
in their calculations. This is not an unexpected result and a 
temperature different to that e stimated from a simp l e adia - 
batic atmosphere was found bv lKieldsen fc Beddind l|l995h . 
The scaling relations between global asteroseismic param- 
eters obtained by Kepler have been investigated in more 
detail and are presented in another paper (Huber at al., in 
preparation) . 

The verified global asteroseismic parameters deter- 
mined using the methods described in this paper will be used 



in subsequent work to obtain estimates of stell ar structural 
param eters {e.g. using the methods described in lStello et al.l 
l2009bl l and constrain models of stellar structure and evolu- 
tion. The asteroseismic results obtained in this and other 
studies of the observations from Kepler have increased the 
number of stars on which solar-like oscillations have been 
detected and characterised by more than an order of mag- 
nitude and illustrates the great achievements made by the 
Kepler Mission and the Kepler Asteroseismic Science Con- 
sortium. 
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APPENDIX A: TABLES OF FITTED POWER 
LAW PARAMETERS 
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Table Al. Values of the i^max - Ai/ (o) and i^max - Amax (b) power law exponents calculated over the ranges in z^max and Tgg 
indicated below, ai : i/max < 750 //Hz; 02 : 750 < fmax < 1500 //Hz; 03 : i^max > 1500 //Hz; 04 : All i^max; 61 : T^s < 5560 K; 62 : 
5560 < Tgff < 6250 K; 63 : T^g > 6250 K. 
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Table A2. Table of s values derived from power law fits and the parameters obtained from a linear fit to s as a function of Tefj. 



